February 27, 2013 18:33 Philosophical Magazine Isotope-efFect-revised 



Philosophical Magazine 

Vol. 00, No. 00, 00 Month 2012, 1-7 



RESEARCH ARTICLE 



(N 

o 

(N 

> 

O 



Boron isotope effect in single crystals of ErNi2B2C superconductor 

Halyna Hodovanets, Sheng Ran, Paul C. Canfield, and Sergey L. Bud'ko * 

Ames Laboratory, US DOE and Department of Physics & Astronomy, Iowa State 
University, Ames, I A 50011, USA 

(^Received 11 October 2012; final version received 00 Month 200x) 



o 
o 



The influence of local moment magnetism on the boron isotope effect of Ta was studied on 
single crystals of ErNi2B2C. Values of the partial isotope effect exponent of aB=0.10it0.02 
and as=0.10it0.04 were obtained based on two different criteria applied to extract T^. No 
significant change in the partial isotope effect exponent compared to the ones obtained for 
LuNi2B2C was observed. Based on this result we conclude that pair-breaking due to the Er 
local magnetic moment appears to have no detectable influence on boron isotope effect of T^. 



Keywords: isotope effect (IE), borocarbides 



o 

(N 
> 



en 

o 



X 



1. Introduction 

One of the accepted experimental probes of the extent to which superconductivity 
is phonon-mediated is the isotope effect [l|, I4]. In the classical form of the Bardeen- 
Cooper-Schrieffer (BCS) theory, the isotope coefficient a (defined from Tc oc M~°, 
where M is the mass of the element) is equal to 1/2. More realistic theories con- 
sider slight deviation from a = 1/2 \^, iJ]. Measurements of even partial isotope 
effect (for one of the elements in the compound) were instrumental in addressing 
the mechanism of superconductivity in novel superconductors, e.g., non-magnetic 
borocarbides [s, '(?] and MgB2 [7]. 

A significant deviation of the isotope coefficient from the BCS value was theo- 
retically predicted for the case of superconductors containing noninteracting para- 
magnetic ions |8h11||. According to these theoretical views, the introduction of the 
magnetic scattering lowers Tc, and the increase of the pair-breaking scattering rate 
due to the paramagnetic impurities also increases the absolute value of the isotope 
effect exponent 



a 
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^ln(Te) ^ 

"dln(M) I- pilj'{p + 1/2) 



(1) 



where ^{x) is the digamma function, tjj'{x) is the first derivative of the 
polygamma function with respect to its full argument, p = l/27rrp(l -|- X)kBTc, 
Tp is the scattering time that magnetic impurities are characterized by and ks is 
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the Boltzmann constant [12!]. Here a^ is the isotope-effect exponent in the absence 
of the pair-breaking paramagnetic impurities. 

Figure 1 (after [13]) schematically shows the normalized isotope effect exponent 
in the case of magnetic pair-breaking as a function of normalized temperature, 
Tc/Tco, where Tc and Tcq are the critical temperatures of the superconductor with 
and without impurities. 

So far, the number of superconductors with paramagnetic impurities in which 
the isotope effect exponent was evaluated experimentally is rather small. In this 
work we report the measurements of the partial, boron, isotope effect exponent in 
single crystals of ErNi2B2C magnetic superconductor. 



Since their discovery [I3l4l5l]. the rare-earth nickel borocarbides, i?Ni2B2C (where 
R=Y, Gd - Lu), have become a toy box for solid state physicists [1^] due to the 
richness of physics one can find within the series, ranging from coexistence of 
magnetism and superconductivity to complex metamagnetism and to a classicall 
heavy fermion behavior IGHlg]. 



-RNi2B2C crystallizes in a structure derived from the ThCr2Si2 structure, space 
group I4/'mmm., that consists of alternating layers of Ni2B2 (in which Ni is tetra- 



hedrally coordinated by four boron atoms) and sheets of RC [19|, |20|]. Phonon- 
mediated BCS superconductivity was proposed for borocarbides and supported, 
e.g., by specific heat measurements [2l|, however the anisotropy of the supercon- 
ducting gap and the details of the superconducting pairing are still under discussion 



22,123]. 



The phonons responsible for the superconductivity in these materials are thought 
to be the high-frequency boron Aig optical modes [2414271]. For non-magnetic, su- 



perconducting members of the series, the boron isotope effect (IE) of Tc has been 
measured: significant boron IE was reported for polycrystalline YNi2B2C, partial 
isotope effect exponent, aB=0.25±0.04, [5|; and for single crystals of YNi2B2C, 
aij=0.21±0.07, and LuNi2B2C, aB=0.11±0.05 [6]. The exponents for YNi2B2C 
and LuNi2B2C are different due to different molar masses, Debye temperatures 
and electron-phonon coupling constants of two compounds (Y has a molar mass of 
88.9 g/mole and Lu has a molar mass of 175 g/mole; YNi2B2C has a molar mass 
of 240 g/mole and LuNi2B2C has a molar mass of 326 g/mole). 

Several i?Ni2B2C compounds {R = Dy - Tm) exhibit coexistence of supercon- 
ductivity and magnetism. Our goal in this work is to explore whether the boron 
IE of Tc is affected by the rare earth magnetism in the borocarbides. 

To check these predictions in the borocarbides family, boron IE measurements 
in lightly Gd-doped LuNi2B2C or YNi2B2C [23] appear to be a suitable choice. 
Realistically though, with a very fast decrease of Tc on Gd-doping, dTc/dx Ri 50 K 
(that translates to ^ 0.1 K decrease in Tc with 0.2% of Gd), very precise control 
of Gd concentration in two batches (one with ^"B, another with ^^B) may be very 
difficult to achieve. 

For this work single crystals of stoichiometric ErNi2 ^''B2C and ErNi2 ^^B2C were 
chosen. The molar masses of Er and Lu are almost the same, therefore we expect 
Debye temperatures and phonon spectra of ErNi2B2C to LuNi2B2C to be similar 
as well. Tc of ErNi2B2C is 10.5 K, and long-range antiferromagnetic order occurs 
below T/v=5.9 K [29]; superconductivity appears in the paramagnetic state well 
above (by almost a factor of two) the antiferromagnetic ordering. In addition, the 
suppression of the superconducting transition temperatures for pure RNi2B2C (R 
= Ho - Tm) appear to follow the de Gennes scaling [16], similar to what is expected 
for paramagnetic impurities. All of these observations point toward ErNi2B2C being 
an excellent test case for examining possible local moment corrections to the IE of 



Tc. 
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Figure 1. Normalized isotope effect exponent a/ao vs Tc/Tco for Abrikosov-Gor'kov (AG) pair-breaking 
(after [12l|'). ao is the isotope-effect exponent in the absence of the pair- breaking. Filled point indicates 
expected value for ErNi2B2C under assumption of the AG pair breaking contribution. 



The ratio of Tc of ErNi2B2C to that of LuNi2B2C is 0.66 that means that iso- 
tope effect exponent is expected to be ~ 1.6 times higher for ErNi2B2C than for 
LuNi2B2C if Eq. (1) is apphed in this case. Given the error bars in ob reported in 
ref. [6(1, it appears to be challenging but still possible to distinguish between these 
two exponents. 



2. Experimental methods 

Single crystals of ErNi2 ^^B2C and ErNi2 ^''B2C were grown using the flux growth 
technique [291 . |30| with the only exception being that the ingots of ErNi2B2C and 
Ni2B (flux) were arc-melted together in a ratio of 1:2.25 respectively (we found 
this ratio to give the best yield of crystals) before being placed in a vertical tube 
furnace. Same boron isotope was used for both ErNi2B2C and Ni2B ingots. To 
ensure that extrinsic effects are not influencing Tc, the same furnace and the same 
temperature profile were used to grow samples containing ^"B and ^^B isotopes. 
The samples grow as plates with the c axis perpendicular to the large flat surfaces 
of the samples. This single crystal growth technique has proven to be robust and 
reliable in the course of over a decade 29|, l32l43 



To ensure sharp superconducting transitions [3l|], the samples were annealed. 
Nine samples of ErNi2 ^^B2C and nine samples of ErNi2 ^'^B2C with the masses in 
a range from 2 to 25 mg were chosen for annealing. Each sample was wrapped into 
the Ta foil, then all samples, for each isotope separately, were placed into a fused 
silica tube that was continuously pumped on by the turbo-molecular pump to the 
pressure of 10~^-10~ Torr. The tube was heated to 950*^ C initially for 100 h. The 
choice of the annealing temperature was based on Ref. [31|]- After the first 100 h, 
the furnace was turned off and the annealed samples were furnace cooled to the 
room temperature. Subsequently, the same samples were annealed for additional 
200 h, following exactly the same procedure as a lOOh annealing, and leading to 
a total annealing time of 300 h. Each sample was assign a label, from a to i, and 
was tracked individually. This annealing yield the residual resistivity ratio values 
(for in-plane direction of the current) RRR=29 for ErNi2 ^^B2C and RRR=30 for 
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Figure 2. Temperature dependence of normalized magnetization data, M(T)/|M (4.6i<')|, of single crystals 
of ErNi2 ^^B2C and ErNi2 ^"620 annealed for 300 hours. Each curve comes from a different sample. 
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Figure 3. Tc with standard deviation for ErNi2 "'^''620 and ErNi2 "'^"'^B2C single crystals annealed for 300 
hours. 



ErNi2 B2C. This improvement in RRR has been attributed to decreasing the 



number and density of dislocation lines in the crystal 35}]. Annealing also made 
the superconducting transition sharper. 

Magnetic measurements were carried out in a Quantum Design Magnetic Prop- 
erty Measurement System (MPMS) SQUID magnetometer with the field applied 
perpendicular to the c-axis. Zero-field cooled low field (25 Oe) magnetization data 
were taken on increasing temperature. The magnetization data were normalized 
by the absolute value of the diamagnetic moment at 4.6 K. The superconducting 
transition temperatures, inferred from magnetization data, were used to calculate 
the boron isotope constant as- The criteria of extracting Tc and calculation of 
boron isotope constant as will be explained in the next section. 
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Tabic 1. Experimental results with standard deviation for ErNi2 B2C and ErNi2 B2C single crystals. 

ErNi2B2C # of samples Tc criteria fc(lOB) (K) fc("B) (K) AT/^ (K) ag 

annealed 300 h 9 with ^^^B, 9 with "B onset in M/H 11.17 ±0.01 11.06 ±0.02 0.11 ±0.03 0.10 ±0.02 

max in ^^^r^ 11.01 ±0.03 10.91 ± 0.03 0.10 ±0.04 0.10 ±0.04 

3. Data and Discussion 

The normalized magnetization data for all measured single crystals of ErNi2 ^^B2C 
and ErNi2 ^'^B2C annealed for 300 h are shown in Fig. 2. There is a small but observ- 
able isotope shift in the superconducting transition temperature AT/^=Tc(^*^B)- 
rc(^^B) between the data for ErNi2 ^^B2C and ErNi2 ^°B2C single crystals. To 
evaluate AT^^ , the following two criteria for T^ determination were utilized: 1) 
onset in M/H and 2) maximum in ^ ^^ ' . These criteria allowed for the deter- 
mination of the superconducting temperatures that are easily comparable among 
the samples. The data, with the standard deviation, for annealed for total of 300 h 
samples of ErNi2 ^^B2C and ErNi2 ^'^B2C, are presented in Table 1. T^, with stan- 
dard deviation, versus annealing time is shown in Fig. 4 for better visualization. 
The inferred values of AT/ , 0.10ib0.03 K and 0.10ib0.04 K, are equal, within un- 
certainty, for both criteria and are smaller than the ones for LuNi2B2C (0.23ib0.12 
K and 0.16±0.07 K) ^. 

To evaluate the possible statistical errors due to the measurement procedure and 
instrument the susceptibility of one sample of ErNi2 ^'^B2C annealed for 300h was 
measured consequently nine times and Tc of the sample was estimated using two 
criteria mentioned above. Tc estimated from: 1) onset in M/H was (11.03ib0.01) 
K and 2) maximum in ^ ^^ ' was (11.15ib0.004) K. This analysis indicates that 
the statistical errors are much smaller than the ones from the sample to sample 
scatter. 

Tc values given in Table 1 were used to calculate partial boron isotope effect 
given by equation (2) 

Aln(r,) 
"^ = -Aln(M^)' ^^^ 

where Aln(Mij) = In(ie) and Aln(re) = ln(^|^) 0. 

The partial boron isotope effect exponents for 300 h annealed samples are 
aB=0.10ib0.02 and 0.10=b0.04 for the two criteria mentioned above respectively. 
The partial boron isotope exponent calculated for the 300 h annealed samples of 
ErNi2B2C is equal, within uncertainty, to that of LuNi2B2C [6[ despite Er being 
an ion with local magnetic moment. 

The fact that Er bears local magnetic moment and introduces magnetic pair 
breaking into the system seems to have no detectable effect on boron isotope co- 
efficient of Tc which is different from the result of the theory that states that 
introduction of paramagnetic impurities increases the absolute value of isotope 
exponent [121]. It has to be mentioned though that the theory was formulated for 
a very low concentration of the magnetic impurities, so that they can be consid- 
ered non-interacting but the local magnetic moments on Er atoms interact via 
RKKY interaction that eventually is a cause of the long range magnetic order. In 
addition, due to the crystal electric field effect the Er local magnetic moments are 
constrained to the ab plane, that might affect superconductivity. More theoreti- 
cal work needs to be done in order to determine how introduction of interaction 
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among the paramagnetic impurities effects the absolute value of the isotope effect 
coefficient of Tg. 



4. Conclusion 

The influence of the pair-breaking local magnetic moment in ErNi2B2C on the 
boron isotope effect was studied on the single crystals of ErNi2 ^^B2C and 
ErNi2 ^°B2C. Estimated Ar/^=(0.11±0.03) K and (0.10±0.04) K lead to the par- 
tial boron isotope constant aB=0.10ib0.02 and 0.10ib0.04 respectively. The boron 
isotope constants obtained in this work are equal, within uncertainties, to that of 
LuNi2B2C [S] which suggests that the pair-breaking, despite lowering the super- 
conducting temperature, appears to have no detectable effect on the boron isotope 
effect of Tc in this system. 
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